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We analyze the excitation spectrum of a three-dimensionaI(3D) Bose-Fermi mixture with tun- 
able resonant interaction parameters and high hyperfine spin multiplets. We focus on a 3-particle 
vertex describing fermionic and bosonic atoms which can scatter to create fermionic molecules or 
disassociate. For a single molecular level we argue that the low lying excitations of the mean-field 
theory are described by the Fermi liquid picture with a quasiparticle weight and charge which is 
justified by a expansion, expected to be exact in the limit of infinite degeneracy (or very high 

fermionic spin) A^^ — ^ oo. Although we restrict our analysis to one molecular level, our picture is 
expected to apply to dispersive molecular bands. Our emphasis is placed on the novel conditions 
for chemical equilibrium and how many-body chemical reactions renormalize the bosonic chemical 
potential, modifying condensation and superfiuid-insulator transitions. 
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I. INTRODUCTION 

Recent experimental studies of ultra cold atomic gases 
cooled to quantum degeneracy and trapped in optical 
lattice potentials have demonstrated precise control over 
tvifo-body interactions^. The Feshbach resonance, which 
can be exploited by ramping a magnetic field Ai? near 
the resonance Bq, can tune the scattering length ag over 
a wide range of values, generating both repulsive and at- 
tractive interactions (of strength ga,i3 ) between pairs of 
atoms a and j3. The interaction coupling has the follow- 
ing dependence on key quantities, 
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where Saa^ is a field dependent quantity that depends 
on the mismatch of the magnetic moments and a^g is the 
background scattering length. 

Furthermore, optical lattice systems have provided a 
remarkably ideal testing ground for new strongly cor- 
related phenomena. In particular, Bose-Fermi mixtures 
possess a unique physical realization beyond the scope 
of traditional solid state materials. Several experimen- 
tal groups have successfully synthesized mixtures of ^Li 




FIG. 1: Vertices associated with the interaction term Hf,^ 
(a) Disassociation of a fermionic molecule into bosonic and 
fermionic atoms (b) Molecular formation through tunable res- 
onant scattering. 



Fermions with both ^Li Bosons^'"^ and ^"^Na Bosons'*, and 
more recently '*°K Fermions with ®^Rb Bosons^. The 
electrically polar nature of KRb and other heteronuclear 
diatomic molecules allow for novel quantum states driven 
by electric dipole-dipole interactions^, however for sim- 
plicity such effects will not be discussed here. Most recent 
experiments have reported interspecies Feshbach reso- 
nances in mixtures of Rb and Cs atomai. 

Earlier theoretical studies have focused on density- 
density type interspecies interactions. Interesting phe- 
nomena utilizing various mean-phase studies^ has found 
the precise limits of phase separation^rJi. Certain one- 
dimensional approacheai2iH have predicted novel ID or- 
dered phases where strong coupling theories^^ seem to 
break down. More recent numerical work by Lewenstein 
et alJ^ identified several pairing instabilities with a rich 
phase diagram consisting of various paired phases be- 
tween the bosons and fermions, with order ranging from 
CDW to superconductivity. Later studies introduced a 
3-particle interaction which explicitly allows for the cre- 
ation of a boundstate molecular particlei^"— , revealing 
Bose-Einstein condensate depletion and other manifes- 
tations of phase separatio n*^'^" . In this article we shall 
focus on the consequences of this type of interaction. 

Unlike simple electrons, atoms and molecules have in- 
ternal degrees of freedom that build composite particles 
with total spin given hy F = I + S, where / is the nuclear 
spin and S is the electron spin. In the case of fermions, 
the system can exhibit high degeneracy due to hyperfine 
spin multiplets such as the recently observed '^^^Yb Fermi 
gas, with constituents having spin of 5/2^1. Therfore, 
here we consider such a system, consisting of a degener- 
ate Bose-Fermi mixture containing high spin fermions of 
degeneracy = 2F+1, which can combine with bosonic 
atoms via resonant scattering. We derive a Fermi-liquid 
picture of the low-lying excitations of a single molec- 
ular level, having an effective quasiparticle weight and 
broadening width, taking in account non-Gaussian fluc- 
tuations at the order 1 /N^i, . The analysis presented here 
fundamentally unifies the condensed matter physics with 
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atomic physics and chemistry, given that the molecular 
levels which are hybridized through a many-body chem- 
ical reaction can be described through the quasi-particle 
concept. Thus, if the chemical equilibrium conditions 
can be violated, the components of the mixture will be 
driven by a steady reaction rate that will deplete one or 
more atomic or molecular species either by net molecular 
formation or disassociation. 



II. THE MODEL 

Here we study a model of a mixture of spin-polarized 
fermionic atoms and bosonic atoms, both of which can 
be treated as spinless particles. The canonical annihila- 
tion and creation operators are defined for the Bosons 
- b,w and Fermions - /,/^, combined to form fermionic 
diatomic molecules iIj, ip'^ . The most general Hamiltonian 
can be decomposed into the following parts: 

n = no + nbb + nl'} + nlp (2.1) 

where, 

Kbb = 5.. E H:-^J>l-M^' (2-3) 

K = y^j:/kS;J^^' (2-4) 

K = ^T.^^^Mi'+%¥ur^' (2.5) 

k,k' .q 

(2.6) 

The fermionic s-wave scattering is usually suppressed 
due to the Pauli exclusion principle and can be omitted 
from the outset. Furthermore, we take the regime of net 
attraction << such that -ff^J can be neglected. 
The kinetic energy part is given by T-Lq , consisting of dis- 
persion relations that are offset by the respective chemi- 
cal potentials ^. If the system is truly at chemical equi- 
librium then the sum of /i's of the reactants should equal 
that of the products. 

In order to obtain controlled and qualitatively simple 
analytical results when integrating over momenta, we 
shall invoke a wide flat-band approximation where the 
two species of Fermions occupy separate kinetic bands 
with separate Fermi surfaces, and the effective masses 
are large enough such that the dispersion relations are 
slow varying in momentum space near the bottom of the 
bands: el ss eg — /u/, et ~ — fu,^. Although the 
fermionic atoms and molecules occupy decoupled Fermi 
surfaces the hybridization effects of the molecules with 
the atoms will manifest through the broadening of the 
molecular Fermi energy. 

The Bosons are sufficiently condensed into a BEC 
phase such that the total number of Bosons is ~ 
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FIG. 2: (a) The hybridized, one particle Green's function, 
{iijjn) expanded in powers of the bosonic mean field, (b) 
The terms leading to the mixed amplitude G^'J (iujn). 

L^blbQ, for a restricted volume with Eq — fco/2mfc — /if,. 
It follows that the two-body and three-particle Bose- 
Fermi couplings g^^ and g^^ can be tuned by a magnetic 
field via Eq. ljl.ip . thus driving the formation of molecu- 
lar bound states or renormalizing the kinetic bands. The 
vertices associated with the type of interaction are 
shown in FiglT] Term Hbi, describes Boson-Boson inter- 
actions, which affect the condensation fraction Njj/Ni, 
at low temperatures. A sharp transition into the BEC 
phase occurs when the chemical potential reaches a crit- 
ical threshold at which fx^ = E^, where is the low- 
est energy that the Bosons can occupy. This transition 
can be understood as the breaking of the U{1) gauge 
symmetry associated with the change of phase in the 
many-body wave function. For free Bosons in a box, 
fi^ = 0. However, in a ultra cold atomic systems, the 
Bosons are usually confined by a superposition of har- 
monic trapping potentials, resulting in the critical value 
Mb — (^/2)(ti^i + W2 + '^s), where uJi is the trapping fre- 
quency along the zth direction. Hence, in physically re- 
alizable ultra cold atomic systems, the condensate frac- 
tion Njj/Nf, and critical value /ig = E^ strongly depend 
on the details of the confining potentials in addition to 
the 2 and 3 body interaction effects. If the condensate 
fraction is large enough, the Bogoliubov approximation- 
can be employed, resulting in a renormalized Boson mass 
nib and dispersion relation which can be redefined in the 
original Hamiltonian. 

III. MEAN FIELD THEORY 

A. Derivation of the quasiparticle Hamiltonian 

Here we investigate a special case of the system de- 
scribed by the Hamiltonian (|2.6p and restrict our anal- 
ysis to the degenerate molecules at the Fermi surface 
that hybridize with the same number of flat-band atomic 
Fermions. We can express the partition function as a 
functional integral, 

Z = y" VbVb^VfVpVrjjVij!'' exp ^- /:(r)rfr j 

(3.1) 
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with the Lagrangian given by 
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where the constant V is proportional to and the de- 
generacy index m is summed to . Note that the degen- 
eracy index which is attached to the operator fkm reflects 
stoichiometry, in that the number of fermions that partic- 
ipates in the hybridization is fixed by the molecular level 
degeneracy. The resulting model described by Ea. (l3.2p 
has a form that at first glance mirrors the slave boson 
mean-field theory of the infinite U Anderson model, for 
which the mean-field theory and corrections have 
been thoroughly studied. However the operator con- 
straint which projects out the bosona^^ is replaced by the 
bosonic repulsion and chemical potential terms. Never- 
theless, the two models exhibit distinct phenomena which 
will be clarified in a later section. 

Consider the effect of integrating out either both or one 
of the / and ip fields. Since the is linear all fields, 
we can make use of the Hubbard-Stratonovich transfor- 
mation. A proper choice of the integration depends on 
dominant low energy scale of each particle species. RefJ^ 
for instance, integrates out only the molecular field "0 as 
follows, 
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where Zq is the non-interacting partition function of the 
fermionic atoms. By approximating this functional in- 
tegral by the maximum value of its integrand for r in- 
dependent bosonic fields, we arrive at our version of the 
mean-field approximation, where the &'s of the low en- 
ergy ground state will be replaced by their expectation 

values: bo ^ (ba^ — r and 6 J = (^o) ~ '^ith r defined 
as a variational mean-field parameter that minimizes the 
free energy T = —ThiZ. Therefore, the following mean 
field Hamiltonian becomes. 



(3.7) 



with the quasiparticle Hamiltonian given by the expres- 
sion, 



k m rn,p 

(3.8) 

with an effective molecule-Fermion hybridization cou- 
pling V = Vr. 



B. Quasiparticle Thermodynamics 



Z^Pj ^ Z^"''Z^^^ J VbVb^VfVpe-^'ff-fo'»-'"'<-^^'i^ 

(3.3) 

where 2^^™* is an irrelevant constant and Z^'-^ is the par- 
tition function of the non-interacting /,6 systems. The 
effective action becomes 



^e//= r^'EA-™^f^+4) ^b'fLdr (3.4) 
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Thus, the resulting attractive 2-body interaction sim- 
ply offsets the repulsive two-body coupling constant g^j. 
This binding energy reaches a threshold at which the 
scattering amplitude vanishes. Beyond this zero energy 
regime are 3- particle processes shown in Fig. [T] which is 
the central focus of this article. 

On the other hand, if the / fields are integrated out 
followed by the if) fields, then the partition function be- 
comes 
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The quasiparticle free energy J-qp can be determined 
from the partition function. 
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where the quasiparticle density of states is given by 
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and the single quasiparticle Green's function Gqp{iuJn) 
can be determined by solving the coupled equations for a 
broadened energy level that often appears in the analysis 
of local magnetic moments in metals. 
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After substituting Eq. ljXTU)) via Eq. ljXTTI) into Eg. (153)) 
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and integrating by parts we have, 
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where the quasiparticle distribution /(w) is given by, 

1 



(3.13) 



and the low temperature behavior follows from a Som- 
merfeld expansion of each Fermi-Dirac distribution in the 
equation above, 



(3.14) 



where H{lo) is a function that vanishes at as w — 0. 
Next we utilize the wide band limit approximation, where 
the integration is taken over a width 2W with a con- 
stant density of states po, and the summation reduces to 
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and the ground state energy {T — 0) evaluates to. 
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with A = A is determined through a minimization that 
solves the equation, 

9E,s (4)' +^' ] f^b , 2.9,,A _ ^ ^g^^^^ 
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Next, we invoke the Friedel sum rule, which relates 
the quasiparticle occupation to the self-energy correc- 
tions and broadening factor^^^, 

n ^ = (Ar/7r)tan~VA/e^) (3.18) 

Consequently, in order to maintain a finite density of 
levels as oo, then A — > 0. Taking this limit in 

Eq. (13.171) leads to the following zero temperature shift of 
the bosonic chemical potential. 
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with the thermal energy scale of the molecules defined by 
Ep = ksTp which is valid provided that <^W. Evi- 
dently, the BEC will be depleted if this correction shifts 
the bosonic chemical potential above its lowest available 
energy. Moreover, the low temperature behavior can be 
obtained by applying the expansion of Eq. p.l4p in the 
wide band limit to yield. 



5^ib{T) ~ (5/ifa(0) 
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Thus, these corrections determine the large 7V^, low T 
behavior as compared to the characteristic temperature 



n 



However, for a proper sense of the behavior near T^, 



it is necessary to go beyond the mean field theory. 



C. Quasiparticle weight and Fermi liquid theory 

In order to formally establish a local Fermi-liquid pic- 
ture for a broadened molecular level, one must determine 
the quasiparticle weight and establish a 1 to 1 correspon- 
dence between the quasiparticles and original molecular 
levels. The propagator of the broadened molecular Fermi 
level is related to the quasiparticle propagator by the re- 
lation, 
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In the context of the theory of the Fermi liquid, = 
plays the role of the wave function renormalization 
factor or quasi-particle weight z. Furthermore, z can 
be related to several conservation laws that arise from 
constraints on the volume of the Fermi surface. Recent 
studies by Powell et al. have dervied such laws given for 
^3p by applying Luttinger's theorem which states that 
the volume of the respective Fermi surfaces is conserved 
even in the presence of interactions^^. These conservation 
laws take the form, 



z = \{h)\^ = Nb-\m^ 



(3.22) 



Additionally, by applying the relation above one can 
then generate conserved charges. Formally, with any con- 
servation law there is an associated quantum number Q 
such that quasi particle charges can be defined. However, 
Q should vanish if the particle-hole symmetry of the sys- 
tem is strictly upheld. More importantly, the validity 
of these conservation laws can be tested via the measure- 
ment of z and Q, which should be directly accessible from 
spectral line broadening measurements and laser absorp- 
tion that can capture velocity and density distribution 
profiles to a high degree of precision^. 
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(b) 
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FIG. 3: Equivalent contributions to the mean-field free energy 
of order A''^. (a) Ring diagram summation and (b) Hartree 
approximation. 
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FIG. 4: Various contributions to However, only diagram 

(a) contributes at order 1/N-^. 



Eq. p.8|) and an interaction part containing the rest of 
the terms. 



A. Hartree approximation 



D. Diagrammatic approach 

An alternative approach, which will be useful in consid- 
ering non-Gaussian fluctuations is to derive the preced- 
ing results diagrammatically from the standard Green's 
function technique^^, involving an expansion in powers of 
the hybridization V . The resulting thermal, one-particle 
propagator with Lorentzian broadening can be computed 
from the series in FiglJJa). Additionally, another prop- 
agator of interest describes the mixed molecule-Fermion 

amplitude {Tf^ pp}^ = G^\^{iujn), which can be de- 
termined by Fig. El^b) as. 
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Next we can proceed to determine the change in mean- 
field free energy of order due to H^^. The mean-field 
equations follow from a perturbation series in the usual 
way involving Feynman diagrams. The first approach 
follows from the linked cluster expansion where the mean- 
field free energy is computed from the sum of the ring 
diagrams as shown in FiglSlJa). The second approach 
follows from the Hartree tadpole diagram as shown in 
Fig. [Sjb), which is made to cancel the external field 
terms, reducing to the desired mean-field equation. In 
both cases we recover Eq. (j3.17|) . 



IV. BEYOND MEAN-FIELD THEORY 

Previously in the mean field solution, h and 6^ were 
replaced by their expectation values. Beyond mean field 

theory we make use of the transformation h = (h^ = 

6 + r and ~ {^0 ~ ^ ' ' ^ticre b and 6^ are static 
operators that obey the Bose commutation rules. After 
making these replacements, the Hamiltonian separates 
into two parts consisting of the mean-field part given by 



The Hartree approximation can be obtained by ana- 
lyzing the tadpole diagram of FiglSl^b), which contains 
a closed loop of the mixed amplitude G^^J{iuJn), repre- 
sented by Eq. p.23p and a static Bose interaction line. 
Here the mean-field equation can be derived by choosing 
the external field r such that it cancels the the tadpole di- 
agram. This condition is satisfied by adjusting the terms 
that are linear in the Bose field to zero as follows: 



fj.br + 2gbbr 
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After carrying out the frequency summation and k inte- 
gration for a wide flat band, we exactly recover Eq. p.l7p 
and the mean-field result for large N^. 



B. Random Phase Approximation 

The available expansion parameter in our model is 
and in order to justify the Gaussian saddlepoint 
(mean-field theory) of the previous section it is neces- 
sary to examine fiuctuations systematically in a l/A^,/, 
expansion. First we can examine several second order 
diagrammatic contributions to the Boson propagator. 
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where Sb is the exact self energy, which can be expanded 
perturbatively in the usual manner by including various 
virtual processes generated by Wick's theorem. Several 
of these self energy contributions are shown in Fig|4]how- 
ever only diagram (a) will contribute at order 1/7V^. A 
proper renormalization of the perturbation series can be 
accomplished by summing a particular subset of the di- 
agrams to infinite order as shown in Fig. [5l This par- 
tial diagram summation is known as the Random phase 
approximation^^, where in our case the bosonic inter- 
action is effectively renormalized by a pseudo-dielectric 
response of the Fermion-molecule polarization bubble, 
which sums to: 



^Frpa - ^ [~*^" + 4 + 



(4.3) 
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FIG. 5: RPA-like corrections to the mean-field free energy 



which has been carefully evaluated by Read(1985) in the 
context of the the slave Boson representation of the in- 
finite U Anderson model^^"— , the details of which are 
shown explicitly in the appendix \^ After taking the 
derivative of Eq. ljASP with respect to A, and adding 
the mean field contribution of Eq. p.l9|) , the shift of the 
bosonic chemical potential to the leading order in l/N.ff,, 
at T = becomes: 



and the irreducible Fermion-molecule bubble of Fig. [5] 
follows the usual frequency summation rules, 

(4.4) 
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SfibiT = 0) = 
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in which the integral of the RHS of Eq. (|4.5p renormal- 
izes the chemical potential. Let us term this integral / 
and focus on its bandwidth W dependence. / can be 
conveniently re-expressed as, 

^ ^ K^nTj r'^- _x dx 

(4.6) 

where k = (iVA)/(7r4) = (iVA)/(7rT^), for yfcs = 1. 
This expression can be integrated by parts, retaining 
the dominant leading order contribution, and dropping 
higher order logarithmic corrections, 

^ W/e% + + K ln(l + M^/4) 

This yields the final expression, 

5^n,{T = 0) ~ 

iV^A^^4 /x^TT W^-4ln(l + Ty/4) 

n wle% - e\le% + Kln(l -f Wle%) 

(4.8) 

For the appropriate low temperature dependence one 
can again apply the Sommerfeld expansion as indicated 



earlier. Consequently, the convergence of the non- 
Gaussian fluctuations is verified and such corrections are 
sufficiently finite that mean field theory is stable at order 
1 jN.^ and the quasiparticle description holds in the limit 
— j> cx). 



V. DISCUSSION: CHEMICAL EQUILIBRIUM 
AND THE SUPERFLUID-INSULATOR PHASE 
TRANSITION 

The existence of fermions will effect the superfluid- 
insulator transition of the Boson atoms at large ghh ■ For 
^2p type interactions this effect was explored by Refael 
and Demler, who found that the superfluid phase was 
suppressed by fermionic screening fluctuations^^. On the 
other hand, Tewari et. al. found a suppression of the 
Mott insulating lobes arising from a bosonic screening 
effecljS^. Consequently, the total effect is complex and 
the phase boundaries can be enhanced in either direc- 
tion. 

For our particular system with finite g^p, the renor- 
malization of /if, has important consequences in different 
physical phenomena which are experimentally measur- 
able. In the preceding analysis we identified a character- 
istic crossover temperature that directly depends on 
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the magnitude of Sp. Thus, by adjusting the molecular 
Fermi level e^, one sets the crossover scale, given that 
our analysis is valid at temperatures T << T^. In ac- 
cordance with the Kondo-like crossover, we expect fit to 
develop a local peak indicating that our approximation 
is breaking down, separating the low and high tempera- 
ture regimes. This enhancement will affect the fi's and 
reaction rates of the other species, given that chemical 
equilibrium implies that the /i's of the components obey 
a rigid sum rule. Since the bosons contain a self interac- 
tion term, at very large repulsion the system is expected 
to undergo a transition into a Mott insulating state, and 
the mean-field phase diagram contains various insulat- 
ing lobes separated by fixed values of the ratio fii,. The 
Mott insulating states are incompressible implying that 
= and the phase boundary is governed by a zeros 
of some function (j){iii,/U). Therefore we predict qualita- 
tively that by adjusting e|, to smaller values, the phase 
boundaries of the lobes will be shifted downwards via the 
enhancement of fib ■ 

As mentioned earlier we have studied a model of a 
many-body chemical reaction. Chemical equilibrium re- 
quires that E^=reac t^n ts 1^^ = Ej=^™d„cts Mj ■ « the 

binding and disassociation is classical in nature, or ■0 ^ 
f + b then conventional chemical thermodynamics in liq- 
uid solution using the law of mass action, relating the 
relative concentrations and reaction rate IC is appro- 
priate, 

^ = M/ ^ e^'^/'^"^ (5.1) 

which is a thermally activated process that limits IC by 
the activation energy However, in ultra-cold atomic 
mixtures with tunable resonant interactions, the sweep- 
ing of a magnetic field near the Feshbach resonance will 
be governed by Fermi's golden rule instead of the law of 
mass action, controlling the bare coupling g^p and sub- 
sequently JC. Hence, this renormalization of the reaction 
rate is clearly a quantum many body problem beyond the 
scope of ordinary chemistry. 

For an accurate physical interpretation, it is important 
to make note of certain differences between the Bose- 
Fermi mixture that we explored here and the slave Boson 
theory of the infinite U Anderson model. Evidently, the 
final expression for the bosonic chemical potential cor- 
rections given by Ea. (|4.8|) does not contain any of the 
problematic infrared divergences of the Kondo problem 
which are related to the subtle restoration of gauge in- 
variance. There are two primary reasons for this. First it 
must be made clear that we do not enforce any Lagrange 
multipliers onto the Hamiltonian (13.81) , as the Bosons and 
Fermions are real physical particles that do not represent 
vacancies and occupied states. A Lagrange multiplier 
that reduces the span of the extended Hilbert space in- 
troduces an additional mean field parameter which renor- 
malizes e^. Therefore, taking an additional derivative 
of the RPA correction of integral (jA5[) with respect to 



the renormalized Sp yields an infrared divergent inte- 
gral. Another difference is that in our unconstrained sys- 
tem, the Bosons or the molecules do not arise from the 
decomposition of the Hubbard operators, Xm^ = ^Am^i 

Xa^m = '0m^^ A:o,o = b'fb, Xm,Tn = V'Jn^m- Such a de- 
composition for a correlated electron systems retains a 
local gauge symmetry b 6e*^ ,b^ we^^ with a simul- 
taneous change i/; — >■ -06*^ ,■0^ — > V^'^e*^. This symmetry 
is absent in our system, therefore there are no local gauge 
symmetry breaking issues. 

One possible deficiency of our unconstrained mean- 
field theory was that we did not enforce a conservation of 
the volumes of the two Fermi surfaces. In particular, the 
conservation laws generated by Powell et al. derives di- 
rectly from Luttinger's theorem and a proper mean-field 
theory that reflects chemical equilibrium should enforce 
these constraints as Lagrange multipliers onto the Hamil- 
tonian, thus generating a more complicated quasiparticle 
theory as observed in the Kondo problem, possibly gen- 
erating an infrared divergence at order 1/N. 

Another important deficiency of our model is from re- 
stricting the interactions to the hybridization of one de- 
generate molecular level with a wide fiat atomic band 
and the excitation spectrum of the Bose- Fermi mixture is 
qualitatively described by a Lorentzian broadened quasi- 
particle with width A. The final result of Eq. (|4.8l) de- 
scribes a bosonic chemical potential with a renormaliza- 
tion controlled by the ratio e^/W, valid in the regime 

<< W, which shall ultimately impact the conden- 
sation behavior. Although we have specifically chosen 
an A'^-fold degenerate at the molecular Fermi surface, 
the analysis can be generalized to any energy or range of 
energies that are likely to hybridize with the atoms via 

In conclusion, we have explored a wide variety of phe- 
nomena associated with the mean-field description of a 
tunable many-body chemical reaction that can be de- 
scribed in terms of a local Fermi liquid picture. These 
features can be probed directly by velocity and density 
distribution measurements via laser beam absorptioni. 
Furthermore, we have demonstrated a unique relation- 
ship between the condition for chemical equilibrium and 
the criteria for incompressibility observed in the Mott in- 
sulating phases of the supcrfluid-insulator transition for 
the Bosons. Although we only considered a single molec- 
ular level it is conjectured that the certain qualitative 
features of the excitation spectrum of Hamiltonian (12.61) , 
still admit a quasiparticle description, as the restoration 
of the momentum dependent self-energy corrections do 
not alter the Lorentzian broadened form of the mean-field 
Green's function given by Eg. p. lip and subsequently the 
effective bosonic action should be renormalized in a simi- 
lar fashion. However, it is necessary to restore the explicit 
momentum dependence in the interactions to study the 
pairing bubbles of the self-energy that may detect screen- 
ing/dielectric effects, and the novel possibility of nesting 
behavior and other Fermi surface instabilities that may 
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-w 
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i?e[e] 



W 



FIG. 6; The contour C of Ea. (|A2[l is deformed to the contour 
C', allowing for a precise evaluation. 



indicate a breakdown of this quasiparticle picture. 
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Refs,2^"^. Starting with Ea. (|44)) . where Wi. = i^v + 
1)1^1^, ujn = 2TTn/l3 for integer values of v and n, the 
integration can be expressed as a contour integral by 
applying the Matsubara summation procedure. We can 
approximate the atomic Fermion propagator as, 



'i7rpsgn(a;„) 



(Al) 



which for a constant density of states of width 2W yields 
the following contour integral, 



0{W - €)9iW + e)sgn(Im[e] ^ 
e - + iAsgn(Im[e]) 



-de 



(A2) 

where the the contour C surrounds the poles on the imag- 
inary axis a shown in Fig. |6l The finite bandwidth intro- 
duces cuts that span 2W along the real axis, separated 
by the distance iuJn- The integral can be evaluated ana- 
lytically at T = by distorting the contour C — >■ C which 
surrounds the cuts as shown in Fig. El This evaluates to. 



■In 



Sp + zAsgn(w„ 
p + iAsgn(a;„) 



N^A^ /(4) +A2 



(A3) 



27r 



■In 



1^2 



Appendix A: Evaluation of E['^'(za;„ 



Here we carry out the evaluation of the of the sec- 

(2) 

ond order self-energy term Ej^ (iujn) as discussed in 



Next to get the RPA free energy correction, the above ex- 
pression is substituted into Eq. ()4.3|) , and the Matsubara 
summation over the the Bose frequencies uJn can be per- 
formed by extracting the imaginary part and deforming 
the contour to surround the negative real axis to yield. 



AF, 



RPA 



tan 



w 



(jV^A/7r)(tan-i(A/(4 - e)) - tan-i(A/4)) 
. - (4 + {N.^AM ln{[(4 - ef + A^]'^'/W})_ 



de 



(A4) 



To the leading order in this is simplified to yield. 



AFrpa = — 



de — ( 



e-el 



N,,A 



■In 



e — el 



W 



(A5) 
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